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"Dark Matter" in Accretion Disks 
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ABSTRACT 

Using Spitzer Space Telescope photometric observations of the eclipsing, in- 
teracting binary WZ Sge, we have discovered that the accretion disk is far more 
complex than previously believed. Our 4.5 and 8 micron time series observations 
reveal that the well known gaseous accretion disk is surrounded by an asym- 
metric disk of dusty material with a radius approximately 15 times larger than 
the gaseous disk. This dust ring contains only a small amount of mass and 
is completely invisible at optical and near-IR wavelengths, hence consisting of 
"dark matter". We have produced a model dust ring using 1 micron spherical 
particles with a density of 3 g/ cm 3 and with a temperature profile ranging from 
700-1500K. Our discovery about the accretion disk structure and the presence 
of a larger, outer dust ring have great relevance for accretion disks in general, 
including those in other interacting binary systems, pre-main sequence stars, and 
active galaxies. 

Subject headings: stars:dwarf novae - accretion disks - quasras: general 



1. Introduction 

Cataclysmic variables (CVs) are semi-detached binaries in which a white dwarf primary 
is accreting material from its close neighbor, a low-mass Roche-lobe filling object. The 
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donor stars are often taken to be of the main sequence variety, but recent observational 
evidence frequently suggests that they may be more evolved and/or brown dwarf-like stars 
(e.g., Howell et al., 2006). Accretion disks that form in CVs have been well studied and 
modeled using observations obtained from X-rays to the near-IR spectral regions (Warner 
1995). The derived picture of the "standard disk" model is a geometrically thin disk of 
gaseous material orbiting the white dwarf. Disk models, bolstered by observation, follow a 
temperature distribution, T(r), given by 

T = T*{R/R^\ 

where T* is proportional to the white dwarf mass and radius and the mass accretion rate 
(see Frank et al., 1992). Typical accretion disks produce a spectrum that in gross detail 
can be approximated by a sum of blackbody rings, each of increasing radius, decreasing 
temperature, and increasing emitting volume. The observed spectrum thus has a Rayleigh- 
Jeans-like tail that follows a A -2 / 3 distribution. In specific detail, accretion disks act more like 
plane parallel stellar atmospheres, providing an underlying continuum source and emission 
lines from their "chromosphere". The highly broadened emission lines (due to the rapidly 
rotating disk) of H and He provide the classic spectral signature of the accretion disk in CVs. 
The history of accretion disk observations and models in interacting binaries is discussed in 
detail in Warner (1995). 

During the Spitzer Space Telescope cycle 3 round of observations, we obtained the 
first and, to date, only mid-IR time series observations of an interacting binary containing 
a white dwarf primary and a low mass secondary star. This data set is one of only a 
relatively few Spitzer time series observations for any kind of astronomical object. Our 
target, the famous eclipsing interacting binary system WZ Sge, was chosen as a follow on 
to our cycle 2 Spitzer observations of magnetic CVs (polars). Our Spitzer observations of 
magnetic interacting binaries, in which we discovered the presence of extended circumbinary 
dust disks are discussed in Howell et al. (2006), Brinkworth et al. (2007), Hoard et al. 
(2007). WZ Sge, an eclipsing, non- magnetic, interacting binary, was observed in order to 
see if circumbinary dust disks were common in all types of interacting binaries or only those 
with strongly magnetic white dwarf primaries. 

Our Spitzer observations are detailed in the next section. We then present evidence 
using nearly simultaneous optical and near-IR ground-based observations to show that WZ 
Sge was at quiescence and behaving normally at the time of our Spitzer mid-IR observations. 
Next we take a closer look at the phased mid-IR light curves we obtained and use these to 
develop a model for the intra-binary dust torus we believe to be present. 
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2. Observations 

Table 1 provides an observing log for all of the WZ Sge observations discussed in this 
paper. 

2.1. Spitzer Mid-IR Observations 

WZ Sge was observed simultaneously at 4.5 (channel 2) and 8.0 (channel 4) microns 
using the Infrared Array Camera (IRAC; Fazio et al. 2004) instrument on the Spitzer Space 
Telescope (Werner et al. 2004). The observations took place on 2007 July 03 UT (HJD 
2454284.891 to 2454284.958) and lasted for 90 minutes, just over one orbital period for 
WZ Sge. We obtained continuous coverage of WZ Sge using 12 second exposures in both 
channels. As recommended by the Spitzer Science Center for time series observations, we 
did not perform dithers or other offsets during our exposures, achieving the desired time 
series of multiple exposures using "repeats". The IRAC pmasks for channels 2 & 4 showed 
no bad pixels at the center of the array, so we centered WZ Sge in the field-of-view. The 
IRAC images were reduced and flux-calibrated with the S12 version of the Spitzer IRAC 
pipeline and downloaded from the archive as basic calibrated data (BCD). The BCD images 
were corrected for array location dependency using the correction frames provided by the 
Spitzer Science Center. Image mid-times were generated by converting spacecraft modified 
julian date (MJD) from the BCD image headers to heliocentric julian date (HJD). This 
time conversion utilized the algorithm that will generate HJD times in all IRAC BCD image 
headers starting with the S18 processing pipeline. 

Taking these processed BCD images, we then performed IRAF0 aperture photometry on 
WZ Sge and four nearby stars to use as comparisons. These four nearby stars were also those 
used to photometrically calibrate WZ Sge in the optical (see §2.2). We used a star aperture 
size of 3 pixels for WZ Sge as well as all the comparison stars. From experience we found 
that the point response function (PRF) for the IRAC arrays is not well modeled, and we 
obtain more consistent results using aperture photometry rather than via PRF fitting. This 
is also the current advice from the IRAC Instrument Support Team at the Spitzer Science 
Center. Sky subtraction was accomplished by using a 3-7 pixel sky annulus around each 
star. The appropriate aperture correction from the IRAC Data Handbook was applied to 



1 The Image Reduction and Analysis Facility is maintained and distributed by the National Optical 
Astronomy Observatory. 

2 See http : //ssc . spitzer . caltech . edu/irac/dh/ 
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our measured flux densities. We did not perform a color correction other than to utilize the 
isophotal effective channel wavelengths (see Table 4) during our subsequent interpretation of 
the data. This accounts for all but < 1% of the color correction. In addition, we examined 
light curves extracted for several field stars covering a range of brightness and found that, at 
the brightness of WZ Sge, the "relaxation" effect that can sometimes cause a "ramping-up" 
of measured flux densities in IRAC time series exposures had a negligible effect, so we did 
not apply a correction for it. 

Representative IRAC images of WZ Sge at 4 and 8 microns are shown in Figure 1 and 
our mid-IR light curves of WZ Sge are presented in Figure 2. Phasing the Spitzer light 
curves to the well known 81-minute orbital period (0.0566878460 days) and the Patterson 
(1998) photometric ephemeris (photometric phase 0.0 = HJD 2437547.72840), we find that 
at the time of primary eclipse - when the white dwarf and optically emitting accretion disk 
gets eclipsed - we see a highly structured eclipse as well. The primary eclipse in the optical 
and near-IR is considered to be caused by the secondary star hiding flux from the gaseous 
accretion disk with photometric phase zero being the time when the secondary most directly 
hides the hot spot. True phase zero (that is, the time of inferior conjunction of the secondary 
star) is offset from photometric phase zero by -0.046 in orbital phase (Steeghs et al., 2007). 
We will use photometric phase 0.0 throughout this paper. 

In contrast to our previous IRAC observations of magnetic interacting binaries (i.e., 
polars), we did not detect a rising SED from 4.5 to 8 microns but one that is approximately 
consistent with an extended RJ-like tail of the well known gaseous accretion disk. However, 
to our surprise, WZ Sge showed deep, broad eclipses in the 4.5 and 8 micron light curves 
(Fig. 2) indicating that the source of the mid-IR eclipsed light which provides a non-negligible 
fraction of the mid-IR emission is not circumbinary, but within the binary itself. 

2.2. Optical and Near-IR Photometry 

WZ Sge was observed on 26-29 June 2007 UT at the Kitt Peak 0.9-m telescope. CCD 
observations were obtained using the default imager at this telescope, S2KB, and light curves 
were collected in Johnson R band with some data obtained in Johnson B band. All images 
were reduced using IRAF and standard practices with aperture photometry performed for 
WZ Sge and three nearby stars used as comparisons. Light curves were constructed for 
the R band time series observations. Both R and B observations were placed on a standard 
magnitude scale using the listed magnitudes for the comparison stars in Henden & Honeycutt 
(1997) and Henden & Landolt (2001). Our four night combined R band light curve is 
presented in Figure 3 where the typical 1-sigma uncertainty is 0.04 magnitude and the out 
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of eclipse scatter is real. The mean out of eclipse brightness of WZ Sge was B=15.33 and 
R=15.11. 

We note that the R band light curve looks similar to previous light curves of WZ Sge, 
showing a "double-humped" structure outside of eclipse which is attributed to the observer 
viewing the broad side of the hot spot twice per orbit (see Skidmore et al., 2000). The 
primary eclipse is also typical for optical observations of WZ Sge: narrow in time (about 
8 minutes long) and ~0.2 magnitudes deep. The 0.9-m photometric observations occurred 
about 1 week prior to our Spitzer observations and we observed nothing unusual about WZ 
Sge in terms of its accretion state, accretion disk, or eclipse profile. 

We also obtained one orbit of simultaneous R and K s band photometry using Andicam 
on the SMARTS/CTIO 0.9-m telescope. These time series observations were made on 7 
August 2007 UT. The R band light curve looked similar to that shown in Figure 3 while our 
K s band observations required significant binning to provide good S/N thus only a single 
average out of eclipse K s measurement could be obtained. WZ Sge's K s band brightness 
(K s =14.05) was consistent with the K s magnitude listed for this star in 2MASS (Hoard et 
al., 2002). Table 2 provides the near-IR brightness history of WZ Sge. Note that the last 
super-outburst by this star occurred in August /September 2001. 

2.3. Optical Spectroscopy 

Optical spectroscopy for WZ Sge was obtained both before and after our Spitzer ob- 
servations. Spectra were obtained at the Kitt Peak 4-m on 3 June 2007 UT, at the Kitt 
Peak 2.1-m on 26-29 June 2007 UT, and at the CTIO Blanco 4-m on 3-5 July 2007 UT. The 
KPNO 4-m spectra covered 3800-5000A and provided only a few spectra of WZ Sge. These 
data looked typical of those presented herein but did not provide significant orbital coverage 
and are not discussed further. The 2.1-m spectra (4000-8000A, 3 A resolution), covered 4 
sequential nights and thus provided coverage of multiple orbits of the star. The CTIO 4-m 
spectra spanned 4000-7600A at 3A resolution and covered a number of orbits of the star. 
All the spectral observations were reduced in the usual way using standard IRAF spectral 
reduction packages. Observations of calibration frames (bias and flat) were obtained each 
night with arc exposures obtained as needed, near in time and sky position to the target. 
Spectrophotometric standard stars were used to flux calibrate the data. 

A representative mean 2.1-m spectrum is presented in Figure HJ All the Balmer lines 
are in emission and present the typical double peaked structure as usual for WZ Sge in 
quiescence, and are as expected for a high inclination cataclysmic variable containing an 
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accretion disk. Figures [5] and [6] show the orbitally modulated S-wave, that is, the tell-tale 
signature of the gas stream interacting with the accretion disk (the hot spot). We note that 
the S wave of the Ha line is accompanied by the S-wave of the Hel 6678A line. Skidmore et 
al. (2000) obtained higher resolution (0.4A), orbit-resolved optical spectra of WZ Sge and 
also presented trailed spectra of the Ha and H/3 lines. Their spectra were obtained prior to 
WZ Sge's last super-outburst in 2001. The Skidmore et al. results are very similar to those 
presented in our Figures [5] and El A more quantitative comparison of the emission lines is 
given in Table [3] in which we list the FWHM of Ha and H/3, measured in a similar manner 
to the Skidmore et al. (2000) paper. Our FWHM values, for the three nights of Blanco 4-m 
spectra obtained in July 2007, are again similar to those given in Table 6 of Skidmore et al. 
(2000). 

All of our contemporaneous optical photometric and spectroscopic measurements and 
our K s magnitude of WZ Sge are consistent with the star's typical values during quiescence 
as outlined and presented in Skidmore et al. (2000) and references therein. We believe this 
consistency provides compelling evidence which assures us that at the time of our Spitzer 
IRAC observations, WZ Sge was displaying its usual quiescent behavior. 

3. Mid-IR Light Curve Analysis 

At the parallax-determined distance of 43 pc for WZ Sge (Harrison et al., 2004), the 
white dwarf will contribute <1% to the mid-IR flux. The optically observed gaseous accretion 
disk tail, roughly following a A -2 3 distribution for a sum of blackbodies (see Ciardi et al., 
1998), should continue to fall in the absence of other flux sources. Thus, the only contributors 
to the flux in WZ Sge at 4.5 and 8.0 microns that were believed to be present prior to our 
mid-IR observations, were the tail of the gaseous accretion disk and the secondary star. 
Using this assumption, we can formulate a model view of what our Spitzer light curves are 
telling us. 

At mid-IR wavelengths, the primary eclipse we observe must hide a source that is 
greater in flux than the hot gaseous accretion disk or any reasonable secondary star at these 
wavelengths. Additionally, the eclipsed source must be far more extended that the gaseous 
accretion disk observed in the optical due to the total time that the eclipse lasts in the 
mid-IR. The primary eclipse width we observe is nearly 0.3 in orbital phase (~24 minutes), 
compared to the optically observed primary eclipse time of only 0.1 phase (~8 minutes). 

The phasing of the light curve tells us that cool material is being eclipsed by the sec- 
ondary star from phase 0.74 to phase 0.14, with a slight asymmetric distribution to earlier 
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phases. The time scale of the ingress, about 0.15 in phase, suggests a linear dimension of 
the eclipsed material of ~ 120,000 km or about the radius of the white dwarf Roche Lobe. 

The depth of the primary eclipse at 4.5 microns is about three times that observed at 
8.0 microns, suggesting that the eclipsed material emits more flux (by this same ratio) at 4.5 
microns. But, the color difference light curve (Flux [4. 5] /Flux [8.0]), also shown in Figure 2, 
reveals a flat color term throughout the entire orbit. This lack of color variability sets a strict 
condition on the source being eclipsed; it must be generally featureless and monotonically 
decline from 4.5 to 8.0 microns. 

The cool, dusty material being eclipsed must be very optically thin, being composed of 
relatively large (few micron) grains. We do not see such a large eclipse extent in the optical 
or near-IR (1-2.5 microns); thus our Spitzer observations have discovered that the dust disk 
contains dark matter - matter invisible in the optical or near-IR band-passes. The dust ring 
appears to be skewed from a uniform distribution and is concentrated on the side toward 
the secondary star, much like the accretion disk bulge observed in the interacting binary EX 
Hya (Belle et al., 2005). It is likely that the dust ring is not uniform in density, thickness, 
or radius and may show a bias toward the secondary star due to tidal forces (as it is closer 
to the secondary star than the white dwarf along the line of centers or due to the presence 
of the accretion stream and associated impact site that produces the optical bright spot. 

Near phase 0.5, where a secondary eclipse would occui§|, we see a dip in the light curve. 
This could be interpreted as a secondary eclipse caused by the vertically extended accretion 
disk itself hiding portions of the cool dust torus and/or the secondary star. If the cool 
dust ring is not extended, this secondary "eclipse" would be entirely due to the secondary 
star. However, the long duration (0.3 in phase, equal to the primary eclipse length) and the 
asymmetric nature of it, suggest it is not an eclipse at all, but simply a time of lower flux 
from the system components such as the similarly observed optical feature seen by Patterson 
et al., (1998, see the top panel of their Fig. 2). 



3 Note, no secondary eclipse has ever been seen in WZ Sge in the optical or near-IR and none would be 
expected as the binary inclination is only 71 degrees. 
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4. Spectral Energy Distribution Model 
4.1. The Code 

In Brinkworth et al. (2007), we introduced our IR spectral energy distribution (SED) 
modeling code for CVs, which we applied in that work to photometric data for several 
magnetic cataclysmic variables (i.e., polars). Improvements to the code were described in 
Hoard et al. (2007), of which the most significant was the addition of a physical cyclotron 
model component. That component is not needed in our current study of WZ Sge (which 
does not contain a strongly magnetic WD). We have continued to improve the code by 
implementing a model component to represent the accretion disk in non-magnetic CVs. 
This component assumes an optically thick steady state disk composed of concentric rings 
emitting as blackbodies, following the "standard model" prescription in Frank, King, & Raine 
(2002). The accretion disk component (denoted with subscript ACD) is fully parameterized 
by a (constant) height (/i ac d), inner and outer radii (-R ac d,in and -R ac d,out), and a mass transfer 
rate from the secondary star (M). The standard model radial temperature profile from Frank, 
King, & Raine (their Equation 5.41) is used to determine the temperature in each disk ring. 
The inclination of the disk is assumed to be the same as the binary inclination; higher 
inclination systems have smaller projected surface area, resulting in an overall decrease in 
the effective brightness of the disk. This component is discussed in more detail in Hoard et 
al. (2008), where we have successfully applied it to observations of the novalike CV V592 
Cassiopeia spanning 0.1-24 microns. 

The shape of the resulting model accretion disk SED is similar to a slightly flattened 
and stretched blackbody function; that is, it displays a relatively broad peak with a rapid 
drop in brightness on the short wavelength side and a gradual (Rayleigh- Jeans-like) decline 
on the long wavelength side (e.g., see Figure 20 in Frank, King, & Raine, 2002). In general, 
increasing either the disk area (i.e., making -R ac d,out larger) or the mass transfer rate causes 
the overall accretion disk SED to become brighter. Increasing -R ac d,out also tends to shift the 
peak of the accretion disk SED toward longer wavelengths, as more surface area is present in 
the cooler outer regions of the disk. Increasing M tends to shift the peak of the accretion disk 
SED toward shorter wavelengths, as the temperature in the inner disk increases. Thus, an 
accretion disk SED that contributes significantly (peaks) at short wavelengths will contribute 
less and less at longer wavelengths. An accretion disk SED that peaks at long wavelengths 
will contribute little at wavelengths only slightly shortward of the peak. So, in order to 
contribute significantly at long wavelengths (a case of interest here), the accretion disk SED 
must either peak at short wavelengths and be very bright overall, or peak at long wavelengths 
(in which case it will not contribute significantly at short wavelengths). These generalizations 
apply regardless of the exact nature of the accretion disk SED, provided that the disk follows 
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a physically plausible radial temperature profile in which disk temperature decreases as disk 
radius increases. 

We note that numerous attempts have been made to generate more realistic model CV 
accretion disks (e.g., Puebla et al. 2007; Collins et al. 1998); for example, by including effects 
of limb darkening, complex opacity laws, and so on. Typically, these modeling attempts are 
applied to the high M (> 1O~ 1O M yr 1 ) cases (e.g., novalike CVs, dwarf novae in outburst), 
which makes them difficult to generalize to systems such as WZ Sge, which have extremely 
low mass transfer rates. In comparison, our accretion disk representation is an admittedly 
simple model. However, it is intended solely to approximate a realistic continuum shape for 
a steady state optically thick disk. It does not attempt to reproduce, for example, spectral 
line absorption or emission (which would require consideration of optically thin regions, 
such as disk chromospheres or coronae/winds). It has at least been found in more complex 
models (Collins et al. 1998) that the continuum spectrum emitted by an accretion disk is 
independent of the opacity, which provides support for the use of our simple model. Another 
commonality of more complex disk models is the strong degeneracy (i.e., non-uniqueness) 
of the solutions (Puebla et al. 2007). We are already concerned about the number of free 
parameters in our approach to matching the observed SED of WZ Sge (see §4.3). Previous 
detailed models have shown that for the continuum, at least, utilizing the simplest approach 
to representing the accretion disk (and other model components) works well. 



4.2. The Data 

Figure [7] (top) shows the observed SED for WZ Sge, composed of the weighted mean 
average of the out-of-eclipse portions of our IRAC light curve (i.e., excluding orbital phases 
0.8-0.2), the 2MASS photometry, the weighted mean average of the out of eclipse ground- 
based B and R light curves, and a UV point at 1580 A from the Hubble Space Telescope 
spectrum shown in Steeghs et al. (2007). The 2MASS J, H, and Ks magnitudes and the 
KPNO 0.9-m optical photometry (and their uncertainties) were converted to flux density in 
mJy using the absolute calibration zero points given in Cohen et al. (2003) and Cox (2000), 
respectively. These data are listed in Table HI 

Our measured K s band brightness from Andicam on 7 Aug 2007, near in time to our 
Spitzer observations of WZ Sge, is in close agreement with the 2MASS value. Consequently, 
we believe that the 2MASS JHK S photometry, which was obtained in 1998 (i.e., before the 
2001 outburst), is still representative of the current post-outburst near-IR quiescent state of 
WZ Sge. 
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The flux density error bars for the IRAC data reflect the weighted standard deviations 
of the weighted averages plus a 3% allowance for systematic uncertainty in the flux calibra- 
tion. The flux density error bars for the 2MASS and ground-based optical data are the la 
photometric uncertainties propagated through the corresponding magnitude to flux density 
conversions. The flux density uncertainty on the UV point is estimated at 10%. In all cases, 
the "error bars" in the x-direction (wavelength) show the widths of the photometric bands. 

4.3. The Model 

Figure [7] (bottom) shows the SED data for WZ Sge with a representative model con- 
sisting of a WD, secondary star, gaseous accretion disk, and dusty ring around the outer 
accretion disk. The parameters of these model components are discussed below. 

The WD (model component subscript WD) is represented by a blackbody curve with 
T wd = 13, 000 K, which was estimated by Godon et al. (2006) as the likely inter-outburst, 
quiescent WD temperature based on a cooling curve constructed during the decline from the 
2001 outburst. The WD has M wd = O.88M and R wd = 6.5 x 10 8 cm, obtained from the 
radial velocity study of Steeghs et al. (2007). The projected surface area of the WD has been 
reduced by 25%, consistent with obscuration by an accretion disk (see below) that extends 
outward from -R ac d,i n = 1-Rwd & t the system inclination of 77° (Spruit & Rutten 1998). If 
the full surface area of the WD is used, then the WD component alone accounts for all of 
the observed UV-i? flux density, which is inconsistent with spectroscopic observations in the 
optical, which show a significant non-WD contribution (see Figure 4; also Mason et al. 2000 
and Gilliland et al. 1986). 

The secondary star (model component subscript SS) is represented by an empirical 
template for an L5.0 brown dwarf based on 2MASS and Spitzer observations of the L5.0 
stars GJ1001BC, SDSS J053951.99-005902.0, and 2MASS J15074769-1627386 (Patten et al. 
2006). As noted by Steeghs et al. (2007), their mass of M ss = 0.078M o for the secondary 
star in WZ Sge implies a spectral type of ~L2; however, (as noted by those authors) this is 
contradicted by the lack of detection of the secondary star in the near-IR (e.g., we note that 
the .fr s -band flux density of an L2 brown dwarf at d = 43.5 pc by itself exceeds the observed 
K s value for WZ Sge by almost 50%, whereas the K s -b&nd flux density of an L5 brown dwarf 
is only fa 40% of the observed K s value). 

For the accretion disk (model component subscript ACD), we assumed that the disk 
starts at the WD surface -R a cd,in = 1-Rwd- The outer radius of the accretion disk is constrained 
by the relatively sharp peak in the UV-optical bands. If -R ac d,out ^ 3i? w d, then the peak of 
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the accretion disk SED is shifted too far to the red, resulting in a much broader peak in the 
total model SED than is observed. This is a small accretion disk, but not unreasonably so 
considering that WZ Sge went through a major outburst in 2001 which likely completely (or 
nearly so) disrupted the accretion disk. The estimated mass transfer rate in Tremendous 
Outburst Amplitude Dwarf Novae (TOADs) like WZ Sge is M< 10 _11 M Q yr" 1 (Howell et 
al. 1995). At low TOAD mass transfer rates, it could take a long time for the disk to 
re-establish itself. In the case of WZ Sge, the mass transfer rate is constrained to be even 
smaller (~ 10~ 12 M Q yr _1 ) than the TOAD maximum so that the observed data at short 
wavelengths are not exceeded by the combined WD and accretion disk model components. 
The steady state accretion disk model component has a resultant maximum temperature of 
w 9200 K, sufficiently hot to account for a significant amount of ionized H, as required by 
optical spectroscopic observations of WZ Sge. We note that forcing the accretion disk SED 
to account for the IRAC flux density values causes the total model SED to far exceed the 
observed flux densities at shorter wavelengths. 

To account for the remaining long wavelength (IRAC bands) flux density, we utilized 
a circumstellar dust ring surrounding the accretion disk (model component subscript CSD). 
This component is calculated following the prescription in Brinkworth et al. (2007) and 
Hoard et al. (2007); namely, it is an optically and geometrically thin annular disk composed 
of dust grains with r = 1 jum and p = 3 g cm -3 . The dust ring is divided into 1000 concentric 
rings following a T cx R~ 3 ^ profile. Each ring contains the same mass of dust, resulting in 
a decreasing dust density at larger radii. The model is parameterized by picking inner and 
outer radii (which control the range of dust temperature and, hence, the overall shape of 
the SED), a height (which, with the two radii, determines the overall volume of the dust 
ring), and a total dust mass (which determines the number of re-radiating dust grains and, 
hence, the overall brightness of the dust ring). We constrained the inner radius to have a 
value such that the maximum temperature in the dust ring is in the range 1000-2000 K, 
corresponding to the likely temperatures of dust grain sublimation. The inner boundary 
condition for the radial temperature profile is provided by the WD temperature and radius; 
in principle, the accretion disk also heats the dust, but we have not accounted for this, 
primarily because of the complexity of calculating heating of the dust due to two sources, 
as described in Brinkworth et al. (2007). We note, however, that the end result of having a 
higher boundary condition temperature would be to move the dust ring radii to larger values 
in order to produce the same inner and outer dust temperatures. The model SED shown 
in Figure [7] requires inner and outer radii of 11 and 30 i? w d, respectively, resulting in dust 
temperatures of < 1500 K throughout the ring. The requisite total dust mass is ~ 3 x 10 17 g 
(about 1.5 x 1O~ 16 M or about 4 x 10- 9 M moon ) This mass is equivalent to a relatively small 
asteroid in our solar system (e.g., about 5% of the mass of 433 Eros; Baer & Chesley 2008), 
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which suggests a possible origin for the dust in the tidal disruption of an asteroid by the 
WD. Alternatively, the dust could have condensed out of material ejected from the WD and 
accretion disk during the 2001 outburst that did not escape from the WD Roche lobe. 

All of the model components were scaled to a distance of d = 43.5 pc, which was 
determined from Hubble Space Telescope Fine Guidance Sensor trigonometric parallax mea- 
surements (Harrison et al. 2004). The values of the model component parameters are listed 
in Table [5J Figure M shows a to-scale diagram of the system. 

As detailed more fully in §3.2.5 of Hoard et al. (2007), the model SED shown in this 
work is not, strictly speaking, a unique solution, in the sense that very similar results can 
be achieved from somewhat different combinations of input parameters. We have tried to 
minimize this effect as much as possible by constraining plausible parameter ranges based on 
whatever other information, observational data, and reasonable assumptions are available. In 
some cases, we can set firm upper and/or lower limits on the plausible ranges for parameters. 
For example, the inner and outer radii of the circumstellar dust ring are bounded by the 
radii at which the ambient temperature exceeds 2000 K (and dust sublimates) and the WD 
Roche lobe, respectively. 

Other limits can be set by the requirement that the total model SED (or that of an 
individual model component) cannot (greatly) exceed the nominal observed SED - as de- 
scribed above, the mass transfer rate cannot be significantly larger than ~ 1O _12 M yr _1 
or the WD and accretion disk together are brighter than the observed optical-near-IR data. 
Similarly, we can exclude L2 and earlier spectral types for the secondary star in WZ Sge 
because they are too bright at certain wavelengths compared to the observations. We note, 
however, that later spectral types for the secondary star (which, regardless of the exact value 
of its mass, can be largely a matter of the temperature - hence, cooling time - of this low 
mass, brown-dwarf-like object) cannot be strongly constrained. In fact, an L7 secondary 
star would improve the model results because the "hump" in the L5 SED between m 2-3.5 
/im (which is actually unabsorbed continuum on the long wavelength side of a prominent 
methane absorption band) shifts to longer wavelengths in the L7 SED. This makes repro- 
ducing the sharp bend in the near-IR JHK S bands easier. Nonetheless, we have used an L5 
secondary star as the most conservative estimate for WZ Sge. 

On the other hand, the interplay between some model parameters (e.g., the dust grain 
size and density) is such that the effect of changing one parameter can be exactly offset by 
adjusting another. In such cases, we have had no choice but to fix the model parameters at 
reasonable values. In Table 5, we have included an estimate of the range over which certain 
model parameters (e.g., those which are not tightly constrained by published observations 
and analysis) still produce a total model SED that, while not as good as the representative 
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model shown in Figure 7, nonetheless does not deviate significantly from the observed data 
(i.e., it reproduces the approximate flux density levels and gross details of the observed SED; 
namely, a narrow peak at optical wavelengths, a plateau in the near-IR out to 
and a gradual decline at longer wavelengths). We have made an attempt to re-optimize the 
total model SED after changing the value of one parameter by adjusting the others (e.g., 
the effect of decreasing the mass transfer rate can be somewhat offset by increasing the 
accretion disk outer radius). The accretion disk and circumstellar dust ring components 
are largely detached from each other in this regard: each of them is most significant over a 
different wavelength range (short wavelength for the accretion disk, long wavelength for the 
dust ring). This means that the largest discrepancies between model and observation, when 
changing the parameters of either of these components, occurs in the near-IR, where they 
overlap. 



5. Discussion 

Using Spitzer Space Telescope observations, we have discovered evidence that a cool 
circumstellar dust ring surrounds the hot, gaseous accretion disk in WZ Sge. This dust ring 
lies inside the binary within the Roche lobe of the white dwarf. Our best model suggests that 
the dust ring contains ~ 3 x 10 17 g of material, modeled as 1 micron spherical dust grains 
with density of 3 g cm -3 (silicate). Reasonable changes to the total mass of dust and/or 
the grain size and density can be balanced by adjusting the other parameters to produce 
equivalent model spectral energy distributions. The dust torus is optically thin and unseen 
in optical and near-IR bandpasses because of its low temperature and particle density, thus 
consisting of "dark matter" . The dust ring is believed to extend from beyond the outer edge 
of the gaseous accretion disk to 30i? w( j, and to have a temperature profile that declines from 
1460 K to 690 K. 

This discovery was a complete surprise as it seems counter-intuitive to expect dust so 
close to a relatively hot star (i.e., the 13,000 K white dwarf). We do know, however, that 
observations have shown that dust exists in the colliding winds of hot (20,000 K to 50,000 
K) Wolf-Rayet stars only a few tens of AU from their surfaces (Pittand & Dougherty 2006; 
Tuthill & Monnier 2008). The apparently high density of material formed at the shock regions 
allows for efficient shielding of the highly energetic photons that would otherwise destroy the 
dust grains. Perhaps the accretion disk contains a high enough outer disk density and a cool 
enough outer disk temperature to provide a shield against destruction of the dust by UV 
photons. Evidence for such outer accretion disk conditions has been observed for WZ Sge in 
the past in terms of the detection of CO and H 2 emission from the accretion disk (Howell et 
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al. 2004). 

Our discovery of circumbinary dust disks in polars (Howell et al. 2006, Brinkworth et 
al. 2007, Hoard et al. 2007) may suggest that all interacting binaries have such material. 
However, if such a circumbinary dust disk existed in WZ Sge with similar emission properties 
as those we observed in polars of nearly equal orbital period, then it would have dominated 
the mid-IR flux distribution in WZ Sge. One important implication of that scenario is that 
we would not have observed a deep (or possibly any) mid-IR eclipse. The presence of the 
prominent mid-IR eclipse at both 4.5 and 8.0 microns argues that the dust must be inside 
the binary. In addition, owing to the much larger surface area of a circumbinary disk, the 
mid-IR flux distribution in WZ Sge would have been much brighter in comparison to its 
optical and near-IR flux densities (e.g., as observed for EF Eri, see Hoard et al. 2007) if the 
dust in WZ Sge was circumbinary instead of circumstellar. 

The origin of the dust is a subject of debate. One possibility is that accretion disk 
material ejected during outburst or superoutburst can become dense enough in shocked 
regions associated with the outburst to both form dust and provide shielding from destruction 
by UV photons for some time period. Material that did not reach escape velocity would fall 
back in toward the white dwarf where its angular momentum would form it into a disk. WZ 
Sge had its most recent superoutburst in 2001 and if the dust formation is associated with 
such events, reobservation at mid-IR wavelengths over time may show a change (decrease) in 
the dust disk. Dust may also form in the atmosphere of the cool brown dwarf-like secondary 
star and be transferred to the white dwarf via the inner Lagrange point. As it sublimates, 
most of it becomes part of the gaseous accretion disk but any grains that remain may be 
the formation particles of the dust disk. Finally, we know of the discovery by the Spitzer 
Space Telescope of dust disks around isolated white dwarfs (e. g., Garcia-Berro et al. 2007). 
These disks are believed to be the remains of progenitor solar systems and are replenished 
by the infall and breakup of cometary or Kuiper belt objects. This mechanism may be true 
in WZ Sge as well, albeit complicated by the binary nature of the system. 

The gaseous accretion disk will behave like an accretion disk and be replenished from the 
(gaseous) matter stream from the LI point. But the dust ring does not act like an accretion 
disk (i.e., no viscous interactions), it acts more like the rings of Saturn. There certainly 
could be (and probably is) some tenuous "grey zone" between the gaseous accretion disk 
and the dust ring, that is part gas and part dust (and dust precursors). How long might 
such dust last if its only destruction mechanism is to spiral in close to the white dwarf via 
the Poynting-Robertson effect and be destroyed through sublimation? The spiral in time for 
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a (spherical) dust grain is given by the expression (Rybicki & Lightman 1979), 

47r P c2 d 2 

where p is the grain density, L* is the luminosity of the white dwarf, R is the grain radius, 
and r is the radius of the initial orbital location. As the grain spirals in toward the white 
dwarf, it is assumed that the orbit is always circular. For the solar system, a one micron 
dust particle orbiting at 1 AU requires about 1000 years to spiral in and sublimate due to 
Poynting-Robertson drag. In the WZ Sge system, we ignore the partial shadowing of the 
white dwarf by the accretion disk, a fact that in reality would lower the effective value of L* 
and increase the spiral in time; r oc L" 1 . For WZ Sge's dust disk, assuming it! ~ 1 micron 
grains, p = 3 g cm -3 , and r is, on average, 15 Rwd, Poynting-Robertson drag will spiral 
the dust into the gaseous disk causing grains to evaporate, in 0.87 years. The source of the 
dust in the intra-binary dust disk must therefore replenish itself approximately yearly. This 
short time scale would seem to eliminate the possibility that the dust forms only via ejected 
material during superoutburst, the last being 6 years prior to our observations. However, 
the mass transfer rate for WZ Sge's donor star is 10~ 12 M yr _1 , while our dust disk mass 
estimate is 10~ 16 M , only 1/10,000 of the mass loss rate each year. Thus, it appears that 
the dust torus can easily be replenished continuously by a small amount of dusty material 
transferred from the cool brown dwarf-like secondary star. Other than Poynting-Robertson 
drag, the dust ring should be fairly static. 

Further study is required to determine if all cataclysmic variables have dark matter in 
the form of a greatly extended dust disk surrounding the gaseous accretion disk. Certainly 
in cases where the white dwarf is very hot (e.g., 40,000-50,000 K as found in some novalike 
CVs), the temperature everywhere inside its Roche lobe is likely to be formally too high for 
dust to exist. However, even in these cases, shielding by the gaseous accretion disk might 
still allow some dust to survive inside the white dwarf Roche lobe. Temporal study would be 
helpful as well to aid in the determination of the origin of the dust. For any accretion disks 
that do have a similar dust ring surrounding them, measurements attributed solely to the 
flux of the central object at mid-IR wavelengths may need to be reexamined. In the case of 
WZ Sge, attributing our measured mid-IR flux densities to the white dwarf, secondary star, 
and accretion disk would have overestimated their contributions by a factor of 2-3. If all 
accretion disks are complicated by the presence of previously unsuspected "dark matter" in 
the form of cool dust, then this could have important implications for the development of disk 
models, accounting for the full energy budget of accreting systems, and the determination 
of bolometric (integrated) luminosities of the components in systems ranging from pre-main 
sequence stars to active galaxies and quasars. 
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Table 1. WZ Sge Observing Log 



UT Date 


Telescope 


Type 


Bandpass 


2007 June 26-29 


KPNO 0.9-m 


Phot 


B,R 


2007 June 3 


KPNO 4-m 


Spec 


3800-50001 


2007 June 26-29 


KPNO 2.1-m 


Spec 


4000-70001 


2007 July 3 


Spitzer/IRAC 


Phot 


4.5 & 8.0 fxm 


2007 July 3-5 


CTIO 4-m 


Spec 


4000-76001 


2007 August 7 


SMARTS/0. 9-m 


Phot 


R, K s 



Table 2. Near-IR Photometric History of WZ Sge 



UT Date 


J 


H 


K 


Ref 


1996 Sept. 23 


14.2±0.2 


13.8±0.2 


13.3±0.2 


Ciardi et al, (1998) 


1998 Oct. 2 


14.877±0.039 


14.535±0.051 


14.019±0.059 


Hoard et al. (2002) 


2007 Aug. 7 






14.05±0.21 


this paper 



Table 3. FWHM Values for Ha and H/3, July 2007 



Date 


line 


FWHM 




2007-07-02 


Ha 


16.52 






17.88 




H/3 


12.43 






13.59 


2007-07-03 


Ha 


18.69 






17.16 




H/3 


13.98 






13.12 


2007-07-04 


Ha 


17.79 



17.54 
H/3 13.73 
13.20 
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Table 4. Spectral Energy Distribution Data 



Band 


Wavelength 


Flux Density 


Date Obtained 




(microns) 


(mJy) 


(UT) 



uv 


0.1580 


0.32 ±0.03 


10 Jul 2004 


B 


0.44 ±0.05 


o n r+0.28 


27-28 Jun 2007 


R 


0.70 ±0.11 


9 RR+0- 25 
/.DD_ 27 


27-28 Jun 2007 


2MASS-J 


i 9 or+0.125 

±.zoo_ 115 


-i o 1+0.07 

i.«i_ 08 


23 Sep 1998 


2MASS-H 


i.662ig:iS§ 




23 Sep 1998 


2MASS-K S 


9 1 KQ+0.141 
Z.1J»_q 13Q 


1 Rfi+0-09 
i - b8 -0.10 


23 Sep 1998 


IRAC-2 


4.493 ± 0.508 


1.55 ±0.07 


4 Jul 2007 


IRAC-4 


7.872 ± 1.453 


0.86 ±0.07 


4 Jul 2007 



Table 5. Model Component Parameters 



Component 


Parameter 


Value 


"Good" Range 


Source 


System: 


Inclination, i (°) 


77(2) 




Spruit & Rutten (1998) 




Orbital Period, P or b (d) 


0.05668784740(28) 




Skidmore et al. (1997) 




Distance, d (pc) 


43.5(3) 




Harrison et al. (2004) 


WD: 


Temperature, T w( j (K) 


13,000 




Godon ct al. (2006) 




Mass, M wd (M Q ) 


0.88 




Steeghs et al. 


(2007) 




Radius, i? wd (cm) 


6.5 x 10 8 




Steeghs et al. 


(2007) 


SS: 


Spectral Type 


L5.0 


L5 or later 


Steeghs et al. 


(2007) 




Mass, M 2 (M Q ) 


0.078 




Steeghs et al. 


(2007) 


ACD: 


Mass Transfer Rate, M (M Q yr" 1 ) 


2.2 x 10~ 12 


1.5-3.0 x 10~ 12 


this work 






Maximum Temperature, T ac d (K) 


9200 


8300-9900 


this work 






Inner Radius, i? ac d,in (Rwd) 


1.0 




this work 






Outer Radius, i? ac d,out (-Rwd) 


2.5 


2-4 


this work 




CSD: 


Optical Depth Prescription 


thin 




this work 






Temperature Profile Exponent 


0.75 




this work 






Constant Height, h cs d (Rwd) 


0.1 




this work 






Grain Density, p grai n (g cm" 3 ) 


3.0 




this work 






Grain Radius, r gra ; n (/zm) 


1 




this work 






Inner Radius, -R C sd.in (-Rwd) 


11 


7.5-15 


this work 






Outer Radius, i? csd ,out (-Rwd) 


30 


20-40 


this work 






Inner Temperature, T csc j ;in (K) 


1460 


1150-2000 


this work 






Outer Temperature, T cs( j.out (K) 


690 


550-930 


this work 






Total Mass, M csd (10 17 g) 


2.94 


1.8-7.0 


this work 





-21 - 




Fig. 1. — Spitzer IRAC images of WZ Sge at 4.5 and 8 microns. Note the typical Spitzer 
point spread functions, especially those at 4.5 microns (channel 1), showing the diffraction 
rings. 
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Fig. 2. — Spitzer Space Telescope light curves of the interacting binary WZ Sge at 4.5 (green) 
and 8 (red) microns. The observations are shown at their full resolution of 12 seconds per 
point and the two points shown on the far left give the one sigma errors for each dataset. Note 
the detailed structure in the eclipse light curves, particularly during eclipse. The bottom 
curve (blue), shows the ratio of the 4 /mi to the 8 /mi fluxes revealing the lack of any color 
dependency in the two light curves. 
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Fig. 3. — R band light curve of WZ Sge obtained at Kitt Peak in June 2007. The phased 
light curve contains observations covering four nights and the typical la error per point is 
0.04 mag. The scatter out of eclipse is intrinsic to the star. 




Fig. 4. — Mean out of eclipse spectrum of WZ Sge obtained in June 2007. The major 
emission features are marked. 
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Fig. 5. — Trailed spectra of WZ Sge for the Ha emission line in June 2007 (left) and July 
2007 (right). 
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Fig. 6. — Trailed spectra of WZ Sge for the H/3 emission line in June 2007 (left) and July 
2007 (right). 
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Fig. 7. — Observed spectral energy distribution (top panel) and model (bottom panel) for 
WZ Sge. The photometric data (see §4.2) are shown as filled circles. The bottom panel 
shows a system model (solid line) composed of a WD (dotted line), L5.0 secondary star 
(short dashed line), steady state accretion disk (long dashed line), and circumstellar dust 
ring (dot-dot-dot dash line). See §4.3 for a discussion of the model parameters. 
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Fig. 8. — To-scale diagram of WZ Sge, including the SED model components (based on the 
parameters listed in Tableland discussed in the text). The diagram shows the L5 secondary 
star (horizontal hatched area), WD (small filled circle) and WD Roche lobe (dashed line), 
accretion disk (solid line around the WD marks the outer edge), and inner and outer edges 
of the circumstellar dust ring (dotted line). The inner and outer Lagrange points (plus 
symbols) and system center of mass (cross symbol) are also shown. The system is depicted 
as viewed from "above" (i.e., with the system in the plane of the sky, equivalent to a system 
inclination of 0°). 



